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Abstract. We report the B band optical observations of an old (r ~ 17.5 Myr) radiopulsar PSR B0950+C . 
obtained with the Suprime-Cam at the Subaru telescope. We detected a faint object, B=27™07±0.16. Within 
our astrometrical accuracy it coincides with the radio position of the pulsar and with the object detected earlier 



by Pavlov et al. (1996) in UV with the ITST/FOC/F130LP. The positional coincidence and spectral properties 
of the object suggest that it is the optical counterpart of PSR B0950+08. Its flux in the B band is two times 
higher than one would expect from the suggested earlier Rayleigh- Jeans interpretation of the only available HST 
observations in the adjacent F130LP band. Based on the B and F130LP photometry of the suggested counterpart 
and on the available X-ray data we argue in favour of nonthermal origin of the broad-band optical spectrum of 
PSR B0950+08, as it is observed for the optical emission of the younger, middle-aged pulsars PSR B0656+14 and 
Geminga. At the same time, the optical efficiency of PSR B0950+08, estimated from its spin-down power and the 
detected optical flux, is by several orders of magnitude higher than for these pulsars, and comparable with that 
for the much younger and more energetic Crab pulsar. We cannot exclude the presence of a compact, ~ 1", faint 
pulsar nebula around PSR B0950+08, elongated perpendicular to the vector of its proper motion, unless it is not 
a projection of a faint extended object on the pulsar position. 
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1. Introduction 

PSR B0950+08 is the fourth in the list of the firs t radio 
pulsars discovered in 1968 (e.g., Bell Burncll 1998). Since 



|Cordova et al. (1989) with the Einstein observatory. 



that time it has been extensively studied in a wide radio 
frequency range, from 0.102 to 10.55 GHz. The pulsar pa- 
rameters (period P and its derivative P, age r, magnetic 
field B, dispersion measure DM , spin-down luminosity E, 
distance d, position a, 6, proper motion /i Q , fig, and paral- 
lax 7r) are shown in Table [TJ So far this is the oldest pulsar 
among ordinary pulsars detected outside the radio range. 

The X-ray counter part of PSR B09 50+08 
was first found by |Seward fc Wang (1988) and 
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* Based on data collected at Subaru Telescope, which is op- 
erated by the National Astronomical Observatory of Japan. 



Later, it w as ob served with the ROSAT by Manning and 
Willmorc ( 1994 ), who analyzed the soft X-ray spectrum 
in (0.08-2.4) keV range. The blackbody (hereafter BB) 
spectral fit gave a temperature Tbb = (2.1 ± 0.6) x 10 6 K 
(90% confidence limits) and a very small radius of the 
emitting region i?BB < 20 m. The pulsar distance was 
assumed to be of about 130 pc based on early radio 
parallax measurements by Gwinn et al. (1986) with the 



VLBI. The power law (hereafter PL, F v oc j/~ q ") spectral 
fit was also acceptable, but, owing to the poor count 
statistics for this faint X-ray object, it resulted in a very 
uncertain spectral index —0.2 < a u < 3.1. Some evidence, 
at 99% significance level, of the pulsed X-ray emission 
with pulsed fraction 0.34 ± 0.18 was found by Wang and 
Halpern (1997) using archival ASCA observations in 
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Table 1. Parameters of PSR B0950+08 (from Taylor et al. 1993, unless specified otherwise). 
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lO" 15 


DM 
cm -3 pc 


0=2000, ^2000 


mas yr _ 


n» 
mas 


r 
Myr 


B E 
G erg s _1 


pc 


253 


0.229 


2.97 


9 h 53 m 09:313(3) d 
7°55'36'.'08(4) 


-2.09(8) 
29.46(7) 


3.82(7) 


17.5 


0.24 x 10 12 5.6 x 10 32 


262(5) 



a The position at the epoch of the Subaru observations, MJD 51930. 
b See Brisken et al. (|2002|) 



Parallax based distance 



d The numbers in parentheses are uncertainties 
referring to the last significant digit quoted, 
e.g., 1.234(56) = 1.234 ± 0.056 



(0.5-5) keV range. The ASCA spectral data were also 
equally well fitted either by the BB with a temperature 
T BB = (5.70 ± 0.63) x 10 6 K, higher than that in the 
ROSAT range, and with Rbb ~ 14 m, or by the PL with 
a v = 0.40 ± 0.22. In accordance with NS cooling theories, 
the temperature of the whole surface of this relatively old 
NS should be less than 10 5 K (e.g., |Yakovlev et al. 1999| ). 
This is much less than observed by both the ROSAT and 
ASCA. The high temperature and small emitting area 
inferred from the BB fit can be explained by the presence 
of small hot polar caps at the NS surface produced 
by the impact of relativistic particles from the pulsar 
magnetosphere, see, e.g., Wang and Halpern (1997) and 
references therein. 

An optical counterpart of PSR B0950+08 was sug- 



gested by Pavlov et al. (1996) based on observations of 
the pulsar field with the ffST/FOC with the long-pass 
F130LP filter (AA = 2310 - 4530 A). A faint (m F i 30 LP = 
27^1) point-like object was found with the 1'.'85 offset 
from the pulsar radio position. The offs et wa s later re- 
vised and decreased to ~ 1" by Pavlov ( 2000 ). If this is 
a pulsar, it is the faintest pulsar ever detected in the op- 
tical. For comparison, the visual magnitude of the young 
Crab pulsar, which is about ten times more distant but 
1.75 x 10 4 times you nger, is 16™65 (e.g., see the review by 
Mignani et al. 2000[). 



Assuming that the d etected optical object is the pul- 
sar, Pavlov et al. ( 1996 ) showed that the extension of the 
ROSAT and ASCA BB fits into the UV-optical range 
gives a flux lower than observed by 3-4 orders of mag- 
nitude. This excludes thermal radiation from the pulsar 
polar caps as a possible source of the optical radiation. 
The assumption that the detected optical flux is due to 
thermal emission from the entire surface of a NS with a 
BB radius Rbb = 13 km yielded the surface temperature 
Tbb - 7 x 10 4 K at d = 130 pc. This is still too high for the 
~ 1.75 x 10 7 yr cooling NS and can be only explained by 
some reheating mechanism operating inside the star (e.g., 
Miralles et al. 1998). On the other hand, the optical flux 



would be in agreement with the extension of the PL X-ray 
fit if the index a u lies within the 0.26-0.35 range. It is in 
agreement with the ROSAT and ASCA data but needs a 
confirmation by deeper observations in X-rays and by the 
detection of the counterpart in other optical bands. 



Possible detection of PSR B0950+08 in the R band 
with the 6m telescope BTA has been reported by 
|Sokolov et al. (1998)| and |Kurt et al. (2000)| . An object 
with R = 25™4 ± 0.3 was marginally (S/N ~ 3) de- 
tected with a short exposure in poor seeing conditions. 
If confirmed, this detection suggests that the pulsar may 
be much brighter in the optical and may have a very steep 
increase of the spectrum towards longer wavelengths than 
one would expect from the detection in the near UV. 

In this paper, we report the observations of the 
PSR B0950+08 field in the B band with the Subaru tele- 
scope. We analyze our data together with the available 
optical-UV data from the HST, and with the X-ray data 
from the ASCA and ROSAT, making use of the recent 
much more accurate radio measurements of the pulsar 
proper motion, parallax, and distance with the VLBA by 



Brisken et al. (2002). Observations and data reduction are 
described in Sec. 2. In Sec. 3 we present the astrometrical 
referencing and photometry, and in Sec. 4 we discuss the 
results and their implications. 



2. Observations and data reduction 

The field of PSR B0950+08 was observed on Jan 21, 
2001 with the wide field camera Suprime-Cam at the 
primary focus of the Subaru telescope. The Suprime- 
Cam ( [Miyazaki et al. 1998j ) is equipped with ten MIT/LL 
2048 x 4096 CCDs arranged in a 5 x 2 pattern to pro- 
vide 34' x 27' FOV with a pixel size on the sky of 
0.201" x 0.201". The pulsar was exposed in one of the 
CCD chips, si006s, with the gain=1.17. Nine 600 s ex- 
posures were taken using the B filter with the bandpass 
close to the Johnson system. The mean seeing was about 
0.7". The 2nd, 7th, and 9th exposures were removed from 



Table 2. Log of observations of the pulsar in the B band. 
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number 


21 Jan 2001 


sec 




arcsec 
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12:41 


600 


1.037 


0.65 


3 


13:21 


600 


1.085 


0.66 


4 


13:54 


600 


1.129 


0.69 


5 


14:08 


600 


1.164 


0.72 


6 


14:22 


600 


1.205 


0.73 


8 


14:50 


600 


1.313 


0.75 
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further consideration because of the problems with the 
telescope guiding system. The observational conditions for 
the rest six exposures are listed in Table 0. 

The standard data reduction, including bias subtrac- 
tion and flatfielding, was performed making use of the 
cedproc task of the IRAF software. To combine the indi- 
vidual dithered images and to get rid of cosmic rays, as 
well as of a trace of an unknown minor object or a dust 
clump with an orbit close to the Earth, crossing our im- 
age E-W approximately 8" away of the pulsar position, 
we made use of the IRAF dither package. A fragment of 
the combined image of the pulsar field with the total in- 
tegrated exposure time of 3600 s is shown in Fig. |I]e 

3. Astrometrical referencing and photometry 

3.1. Astrometry 

The radio position of PSR B0950+08 at the epoch of the 
observations (se e Table was determin ed using the VLA 
observations by Fomalont ct al. (1992) and the most re- 
cent m easur ements of the pulsar proper motion by Brisken 
et al. fl2002| ). 

Astrometric referencing of our image was made using 
positions of five^] of nine reference stars from the USNO 
A-1.0 catalogue visible within the ~ 6' x 13' frame of the 
CCD chip containing the pulsar. The PSFs of these stars 
are not corrupted by the CCD oversaturation effects. We 
used IRAF tasks ccmap/cctran for the astrometric trans- 
formation of the image. RMS errors of the astrometric 
fit are 07086 and 07090 for RA and Dec, respectively, 
whereas the residuals for all stars are < 0713. All of them 
are less than the nominal USNO catalogue accuracy 07 24. 
Combining the RMS errors of the fit, the USNO catalogue 
accuracy and the radio position errors (see Table g) wc es- 
timate total uncertainties of the pulsar position, as well as 
the astrometric referencing accuracy for the whole Subaru 
image, as 0726 in both RA and Dec^. 

3.2. Identification and morphology of the candidate to 
the pulsar counterpart 

As seen from Fig. g, the optical counterpart of 
PSR B0950+08 can only be identified with a faint, 
S/N ~ 7, isolated object clearly visible near the center of 
the 24"x24" fragment of the Subaru image (Fig. 0a). The 
object overlaps with the error circle of the pulsar radio po- 
sition with radius 07 26, corresponding to la uncertainty 
of our astrometrical referencing (Fig. 0c). The source pro- 
file of the object is elongated E-W, with FWHMs - 170 
and ~ 077 along E-W and S-N directions, respectively. 



Image 



available 



FITS format 



at 



http: / / www.ioffe.rssi.ru / astro/NSG/obs / 0950-subaru 
^ U0975.06358239, U0975.06357626, U0975.06357575, 

U0975.06358752, U0975.06358583. 
3 There are no reference stars south of the pulsar within our 

frame and thus we cannot exclude a systematic shift in the 

N-S direction. 



Owing to the low S/N ratio, we do not resolve any point- 
like structure within the source profile. Formal gaussian 
fitting of the whole object profile yields the coordinates 
a 2000 = 9' l 53 m 09!288(17), <5 200 o = 7°55'35"89(26). The 
offset from the radio position is 0/43 W-S. It is within 
1.75c of the astrometrical accuracy and can be considered 
as negligible given our seeing and S/N values. 

To compare our image with the previous optical-UV 
observations of the pulsar field by Pavlov et al. (1996) we 
reanalyzed astrometrical referencing of the HST /FOG im- 
age retrieved from the HST archive. The rotated box in 
Fig. 0a borde rs the FOV of the HST /FOC observations. 
Pavlov et al. ( 1996 ) fou nd th e only point-like object with 
the offset ~1" (Pavlov [2000 ) from the pulsar radio posi- 
tion. Such a large offset makes identification of the FOC 
object with the pulsar and with the Subaru pulsar coun- 
terpart doubtful. We revised the FOC astrometry mak- 
ing use of the FOC position angle and the only reference 
point visible at the north corner of the FOC image, an 
extended object ol (see Fig. 0b). We assumed that ol is a 
distant background object and its proper motion is negli- 
gible. Gaussian center coordinates of ol were determined 
in the Subaru image (a 200 o = 9 ,l 53 m 09!489(2), <5 2000 = 
7°55'4272(1)) and in the FOC image (with the accuracy 
0702), and were used to correct reference frame of the 
FOC image. The coordinates of the FOC pulsar candidate 
in the corrected system at the epoch of the FOC observa- 
tions are a 200 o = 9^53 m 09!298, <5 200 o = 7°55'35776. Given 
that, the discrepancy between the FOC counterpart posi- 
tion and the radio position at the epoch of the HST obser- 
vations, a 200 o = 9 h 53 m 09?314(3), S 20 ao = 7°55'35"88(4), 
marked by "+" in the inset on Fig. 0b, is decreased to 
07 27. This is comparable to the astrometrical referencing 
accuracy of the Subaru image. 

The FOC image is presented in Fig. 0b in the corrected 
coordinates. The overlaid contour map of the Subaru im- 
age shows that the Subaru and HST detected the same 
object near the pulsar position. The FOC contour map in 
Fig. 0c is additionally shifted with respect to the Subaru 
image by —07014 in RA and 07 1 99 in Dec to compen- 
sate for the pulsar proper motion during the 6.75 yr inter- 
val between the HST and Subaru observations. Isophotes 
on the contour maps correspond to the levels (in counts) 
above the background l n = S + na, where S is the mean 
sky value near the pulsar, a is the sky standard deviation 
related to one pixel, n — 1, 2, 3, . . . , 8 and 1, 3, 5, . . . , 15 
for the Subaru and HST maps, respectively. From Fig. 0c 
we see that the Subaru and HST detected the same ob- 
ject near the pulsar position. The better agreement of the 
Subaru and HST source positions after the correction for 
the pulsar proper motion favours this object as the optical 
counterpart of PSR B0950+08. 

We also note that, although the HST/FOC pulsar 
counterpart profile is point-like, the edges of its wings are 
slightly elongated E-W (Fig. 0c). The elongation direc- 
tions in both the FOC and Subaru images coincide and 
are approximately orthogonal to the vector of the proper 
motion of the pulsar marked by the arrow in Fig. 0c. Such 



4 



S.V. Zharikov et al.: Subaru optical observations of the PSR B0950+08 



50 



7:55-f0 



PSR B0950+O8 

SLibiiru Supi irreC?m iiiiT^ 
B 3600s 
Jan 2001 





9:53:10.0 09,6 



36.8 



36.4 



7:55:36.0 



35.6 



35.2 




09.34 9:53:09.30 09.26 09.22 



42 



7:55:40 



38 . 



36. 



34 . 



ft ^ 0<> 





7 9 





PSR 



t* \' V " * o 

1 o _ * l-^Y* 

^ <> 1 

' Pavlov I 
/ HS37FOC 
F430LP 4755s, 
q ♦ Apr1984<l 



9:53:09.7 



09.5 



09.3 



09.1 



Fig. 1. Top Left (^a): 24"x24" fragment of the Subaru 
2001 image of the PSR B0950+08 field in the B band. 
Larger and smaller boxes show the FOV of the HST 1994 
observations (see the right panel) and the borders of the 
bottom panel containing the pulsar, respectively. Right 
7". '4x7'.' '4 HST 1994 image of the pulsar field. The 
pulsar counterpart and an extended object ol, both visi- 
ble in the HST and Subaru images, are marked. The HST 
image is aligned with the Subaru image using ol as a ref- 
erence. The Subaru contour map is overlaid. The pulsar 
close vicinity, ~ 0.6", with the closest Subaru contour and 
the radioposition at the epoch of the FOC observations 
marked by "+", is zoomed in the inset. Bottom ftj^c): 
2"x2" neighbourhood of the pulsar in the B band. The 
HST contour map is overlaid. It is shifted for the pulsar 
proper motion for 6.75 yr passed between the observations. 
The arrow points out the radio position of the pulsar at the 
epoch of the Subaru observations, its length corresponds 
to the pulsar shift due to the proper motion for 10 yr. 
The circle shows la uncertainty of the pulsar position in 
the Subaru image; "+" marks the center of the Gaussian 
fitted to the pulsar counterpart profile in this image. 



an orientation may suggest the association of the elonga- 
tion with a faint torus-like structure of a possible compact 
pulsar nebula as detect ed ar ound young Crab-like pulsars 
(e.g., Weisskopf et al. 2000 ). However, it may be also a 
projection of a faint extended background object at the 
pulsar position. 

To summarize, our analysis shows that the HST 
and Subaru detected the same object. With the al- 
lowance for the pulsar proper motion, the offsets from the 
PSR B0950+08 radio positions at the HST and Subaru 
observational epochs are in the range (0'.'3 - 0'.'4), which 
are within 1.7cr error of the Subaru astrometrical referenc- 
ing and negligible compared to our seeing of Of! 7. 



3.3. Photometry 

Weather conditions were stable during our observations. 
We derived the atmospheric extinction coefficient in the 
B band ks = 0™18 ± 0.02 from the variation of the count 
rates of four stars in the pulsar field with airmass dur- 
ing our observations (see Table ||). Insignificant decrease 
of the extinction, within 1-cr level, was noticed from the 
beginning to the end of observations. 

Photometric referencing was carried out using 
t hree dcfocusc d standard stars from the field PG1047 
( Landolt 1992 ) with V ~ 13, along with five fainter, 
V ~ 21, unsaturated secondary standa rd stars from the 
PSR B0656+14 field flKurt et al. 1998| ) observed in the 
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Fig. 2. Optical and X-ray observations of PSR B0950+08. 
Solid crosses in the X-ray range correspond to the bolo- 
metric fluxes derived from the BB and PL spectral fits 



of the ROSAT (Ma nning and Willmore p94j ) and ASCA 
(Wang and Halpern |l99"7| ) data in the E = (0.08-2.4) and 



(0.5 — 5) keV energy bands, respectively. Square-filled belt 
and brick-filled region show unabsorbed spectra resulted 
from the BB and PL fits of the ASCA data. Stripe-filled 
belt shows unabsorbed BB fit of the HST/FOC data, and 
dashed cross shows corresponding bolometric flux in the 
ROSAT band. The fits are extrapolated toward the X-ray 
and optical ranges, respectively, and their parameters are 
shown in the plot. The widths of belts and region corre- 
spond to lcr uncertainties of the fits. As seen, the optical 
emission is of nonthermal origin and apparently follows 
the extrapolation of the X-ray PL fit of the ASCA data 
with the spectral index a v ~ 0.18 (see Sec. || for details). 

BRI bands the same night. The derived photometric ze- 
ropoint in the B band was 28.28 ±0.02. The instrumental 
magnitudes of the detected optical source were measured 
for a range of aperture radii of (1 — 3) CCD pixels centered 
at the "+" in Fig. 0c. They were corrected for the PSF of 
bright stars (some details on th e method we used can be 
found in Koptsevich et al. 2001). Within the measurement 
errors, the results for different apertures coincided and a 
2 pixel radius was adopted as optimal. The magnitude of 
the detected optical source, B = 27.07 ±0.16, corresponds 
to the absolute flux F B = (5.97 ± 0.88) x 10~ 31 erg cm" 2 
s- 1 Hz" 1 or 0.0597±0.0088 fxJy. The error includes statis- 
tical error of the instrumental magnitude measurements, 
the error of the zeropoint, and an allowance for possible 
atmospheric extinction variations. 

4. Discussion 

In Fig. |^ we combine our results with the avail- 
able data from the optical and X-ray observations of 
PSR B0950+08. In the optical range we neglect the in- 



terstellar extinction, which is expected to be very low for 
this pulsar, E B -v £ 0.02 (e.g., Pavlov et al. |1996|). We 



excluded from our consideration the observations in the 



R band by Kurt et al. (2000) because of very low signal- 
to-noise ratio. In the X-ray range we show unabsorbed 
fluxes resulted from the BB and PL spectral fits of the 



ROSAT (Ma nning and Willmore [99J and ASCA (Wang 
and Halpern 1997) data. The fits are rescaled to the new 
distance value d = 262 ± 5 pc measured by Brisken et al. 



pjQ2D . 

Within the errors the flux Fq coincides with the value 
Ffi3ql p = 0.051 ± 0.003 ^Jy measured by Pavlov et al. 
( |1996| ) with the HST/FOC in the F130LP band. The B 
band (AA - 3944 - 4952 A) considerably overlaps with 
the F130LP band (AA - 2310 - 4530 A). However, their 
pivot wavelengths A B = 4448 A and Ap 130LP = 3750 A 
are different. This difference and close flux values in two 
bands suggest a flat spectrum of the object in a wide 
AA ~ 2300 — 5000 A spectral range. Along with the posi- 
tional coincidence, we consider this flat spectrum, which 
is typical for pulsar optical spectra (e.g., Koptsevich et 
al. 2001; Mignani & Caraveo 200l| ), as an additional ar- 
gument in favour of the detection of the PSR B0950+08 
optical counterpart. 

If the detected object is the pulsar, its flux in the B 
band does not follow the Raylcigh- Jeans law suggested by 
Pavlov et al. ( 1996) ) to explain the optical- UV radiation 
in the only F130LP band as a low temperature thermal 
emission from the entire surface of a reheating/cooling 
NS with Rao = 13 km. In that case in the B band we 
would detect twice smaller flux than measured. This is 
well outside the uncertainties and suggests a nonthermal 
origin of the pulsar emission at least in the B band. 

If the flux in the F130LP band is still dominated by the 
Rayleigh- Jeans tail, we would obtain Tbb = (3.4 ± 0.3) x 
10 5 K, which i s a fa ctor of four-five higher than found by 
Pavlov et al. (|1996|) mainly owing to the change of the 
distance to PSR B0950+08 from 130 pc to 262 pc (see 
the stripe-filled BB belt crossing the F130LP band and 
extended to X-rays in Fig. ||). A NS with such a hot surface 
would produce a bolometric flux ~ (0.8 — 1.9) x 10~ 12 erg 
cm' 2 s _1 in the (0.08 — 2.4) keV band (see a big dashed 



cross in Fig 

io- 14 



3). It is well above the value (2.4 - 7.3) x 



erg cm s , measured by Manning and Willmore 
( |1994D with the ROSAT under the assumption of the BB 
spectrum of the detected X-ray radiation (marked by a big 
cross below the dashed one in Fig. |^) . This means that the 
whole surface of PSR B0950+08 is actually much cooler 
and its emission cannot dominate in the F130LP band. 
Thermal emission from hotter, Tbb = (5.7 ±0.6) x 10 6 K, 
but much smaller polar caps of the pulsar with R pc ~ 
28.6 m, which may explain the detected X-ray radiation 



(Wang and Halpern 1997), can hardly be visible in the 
optical range also because of very small areas of the caps 
inferred from the BB fit of the X-ray data (see square-filled 
BB belt in Fig.|). 

For the above reasons it is most likely that the opti- 
cal radiation of PSR B0950±08 is completely dominated 
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Table 3. Parameters of the radio pulsars and the efficiencies of their optical emission t\b in the B band. 



Source 


logr 


d 


logS 


B 


M B 


log L B 


log??B 




yr 


kpc 


erg s _1 


mag 


mag 


erg s _1 





B0540-69 
Vela 

B0656+14 
Geminga 
B1929+10 
B0950+08 



3.2 
4.1 
5.0 
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T he B band magnitude s are from 



PCTciva^e^^L_^1993) ; 2 Middleditch et al., 1987 (spectroscopic data by Hill et al. (1997)1 give ~ l m .5 smaller magnitude); 



Nasuti et al.. 1997j ; 4 Koptsevich et al. |200l| & |Bignami et al., 1993| ; 6 Pavlov et al. fll996| ) (3a upper limit); 



This work. The distance s d are from 

b |Taylor et al. 1993|, |Finley et alT 



a |Caraveo et al. (2001) 



d Brisken et al. (2002); 



1992 



Anderson et al., (1993)| ; c |Caraveo et al., (1996) 



** The estimates are based on the counterpart detection in the adjacent F130LP band by Pavlov et al. (1996), see Sec. 
Uncertainties of Mb, Lb and t\b include uncertainties of the optical flux and distance measurements. 



by nonthermal emission produced in the magnetosphere 
of the rapidly rotating NS, as it is believed to be for 
young and well studied pulsars like the Crab and PSR 
B0540— 69. Within large uncertainties of the available X- 
ray data for PSR B0950+08 the PL with the spectral in- 
dex ev„ ~ 0.18 matches both the X-ray and optical fluxes 
including the B and F130LP bands. This value is consis- 
tent with what was obtained by Wang and Halpern ( 1997 ) 
from the analysis of the ASCA X-ray data, F PL 



2.77 (+1.33, -1.05)xl0- 13 erg cm" 2 a 
in (0.5 - 5) keV energy range (see the ASCA cross and the 
brick-filled region in Fig. ||). The PL fit of the X-ray data is 
more preferable since the BB fit, which is statistically also 
acceptable, implies by an order of magnitude smaller emit- 
ting area than it is predicted by standard models of hot 
polar caps at the pulsar surface (e.g., Arons 1981 ). The 
inferred spectral index differs from a v ~ 1 which is typi- 
cal for nonthermal soft X-ray radiation of most rotation- 



ASCA — 

0.40±0.22 



powered NSs (e.g., Becker & Triimper 1997). However, it 
can be as low as 0.4 for middle-aged pulsars (Koptsevich 
et al. 2001) and we cannot exclude a decrease of the slope 
of the PSR B0950+08 spectrum towards the optical range 
as it is seen in the case of the Crab pulsar (e.g., Crusius- 
Watzel et al. 2001). More X-ray and optical data are 



needed to check the spectral shape for PSR B0950+08. 
With new data it would be also useful to perform the BB 
+ PL, and/or NS atmosphere + PL fits (e.g., Zavlin et 
al. 1996) to better constrain the parameters of the non- 



thermal and thermal spectral components from the pul- 
sar magnetosphere and polar caps and to estimate their 
contribution to the pulsar emission in different spectral 
bands. 

If the detected object is PSR B0950+08, the pulsar op- 
tical luminosity in the B band, assuming isotropy of the 
emission, is L B = 4:ird 2 F u Aiy = (7.6 ± 1.2) 
s 
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1 , where d,262 = d/(262 pc) is the normalized distance. 
As seen from Table [| and Fig. || (top panel) , where we col- 



Fig. 3. The optical luminosity Lb (top), spin-down power 
E (middle) and optical efficiency in the B band r\B (bot- 
tom) for all radio pulsars detected in the optical range as 
functions of pulsar age r (see also Table ||) . 



lected the data in the B band for all pulsars detected in the 
optical range, Lb of PSR B0950+08 is comparable with 
the luminosity of the middle-aged Geminga pulsar but by 
almost six orders of magnitude smaller than Lb of much 
younger and energetic Crab pulsar. These data indicate 
that, starting approximately from the age of Geminga, 
the decrease of the optical luminosity with the pulsar age 
becomes slower than at early stages of the pulsar evolu- 
tion or even disappears completely, while the spin-down 
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luminosity continues to decrease with approximately the 
same rate (Fig. ^, middle panel). This means that for un- 
known reasons mechanisms of the optical emission may 
start to be more efficient again after ~ 10 6 yr age. The 
ratio of the optical luminosity to the pulsar spin-down lu- 
minosity, r/B = L-q/E, which can be considered as a rough 
characteristic of this efficiency, is plotted in Fig. ^ (bot- 
tom panel). It is seen that optical photons are produced 
by the old PSR B0950+08 with almost the same efficiency 
as by the young Crab pulsar. In this respect it is not sur- 
prising that we see some indications of the presence of a 
compact pulsa r neb ula around PSR B0950+08. Manning 
and Willmore 1994 discuss a compact synchrotron neb- 
ula, which is unresolved with the ROSAT, as the most 
likely explanation of the observed X-ray spectrum of the 
pulsar. The estimated sizes of the possible X-ray nebula 
~ (6" — I') are much larger than ~ 1" we see in our optical 
images. But the difference in the X-ray and optical sizes 
is typical for known compact pulsar nebulae where only 
internal brightest structures of a nebula are visible in the 
optical range (e.g., Weisskopf et al. |2000| ). 

The efficiency derived by Pavlov et al. ( 1 996| ) in the 
F130LP band for another old, - 3.1 x I0 6 yr, PSR 
Bf 929+10 is about (3 - 8) x lQ- 7 d 2 17Q . This pulsar has 
not yet been detected in the adjacent B band, but 
we can assume that its flux in B is close to that in 
F130LP, as it is for PSR B0950+08. Scaling by the fac- 
tor Az/b/Az^fi3olp = 0.24 and by the new radio par- 
allax based distance d = 331 pc (Brisken et al. 2002) 
yields 7 ?B ~ (2.7 - 7.3) x 10" 7 d§ 31 . With this value PSR 
B1929+10 occupies an intermediate position at the rising 
part of the dependence of 773(1"), between Geminga and 
PSR B0950+08, as it is expected from its age. Although 
our consideration of the efficiency evolution is based on 
the data obtained only in the B band, it appears to be 
qualitatively valid for the whole optical range since the 
broad-band spectra of all pulsars are almost flat. Hence, 
the differences of the bolometric optical fluxes of various 
pulsars should not be affected strongly by insignificant 
differences of their spectral slopes in this range. 

Physical reasons for such a high increase of the op- 
tical efficiency at late stages of the pulsar evolution are 
not quite clear. We can only note that the efficiency for 
gamma-ray radiation, rj~ t — Lj/E, also appears to be gen- 
erally higher for older pulsars (e.g., Thompson, 200C ). It is 
difficult to estimate statistical significances of these facts 
and their possible correlation since the numbers of pulsars 
currently detected in the optical and gamma-rays are too 
small (about ten only). Morever, not all known gamma- 
ray pulsars are detected in the optical range (and vice 
versa). However, the observed tendency in both ranges 
seems to be interesting and can hardly be ignored. It is 
obvious that in the both wavelength ranges the radiation 
is nonthermal and originates in magnetospheres of rapidly 
rotating NSs. Thus, the increase of the efficiencies in two 
very different ranges may reflect an overall increase of the 
magnetospheric activity with the NS spin-down. A global 
electrodynamic model of the pulsar magnetosphere with 



the activity caused by the magnetic field-aligned poten- 
tial drop producing electron-positron pairs in the mag- 
netic polar regions of the magnetosphere predicts that the 
efficiency should increase for high energy photons oc P 2 
( ghibata 1995| ). However, the optical data, particularly for 
young pulsars, do not follow this dependence and detailed 
studies of the electrodynamics and radiation processes still 
have to be done to explain the efficiency evolution in dif- 
ferent spectral ranges. 

Further observations of the candidates to the optical 
counterparts of the old pulsars PSR B0950+08 and PSR 
B1929+10 in different spectral bands would be very useful 
to resolve the efficiency problem and to better understand- 
ing the nature of the optical emission of pulsars and the 
evolution of this emission with pulsar age. Measurements 
of their proper motion and the detection of pulsations with 
the pulsar periods in the optical range would be most im- 
portant to provide firm evidence of the pulsar nature of 
the detected optical objects. Simultaneous studies of the 
optical and X-ray pulse profile would provide stronger in- 
dications whether the optical and X-ray emissions are gen- 
erated by the same physical process. Our observations of 
a very faint PSR B0950+08 with the Subaru show that 
new generation of large ground-based telescopes is very 
effective for these studies and could lead to a considerable 
increase of the number of pulsars detected in the optical 
range in the near future. 
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